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Abstract 

We analytically calculate the total and differential cross sections for the pair production of 
doubly heavy baryons in the framework of diquark model. The processes of electron-positron and 
quark-antiquark annihilations are considered. The fractions of doubly heavy baryons in the yields 
of heavy quarks are evaluated numerically. 

1 Introduction 

High luminosities of 5-factories and hadron colliders open a real experimental possibility to observe 
doubly heavy baryons Eqq' and £Iqq>- This prospective stimulates the theoretical interest to study 
the physics of such the baryons: the spectroscopy in the framework of both potential models Q and 
QCD sum rules 0], the lifetimes and inclusive decay modes ||, the mechanism of production in various 
collisions and the rate of yield at accelerators {§. 

The pair production can be essential at the energies close to the threshold, that we consider 
in the present paper. The physical approximation for the calculations are caused by the apparent 
diquark structure of doubly heavy baryons, since the diquark size is essentially less than the radius 
of confinement determining the motion of light quark inside the QQ'q system. We suppose that, at 
first, the calculations of heavy diquark production should be done, and further, the models of diquark 
fragmentation into the baryon || can be applied. 

In paper M the differential and total cross sections for exclusive production of meson pairs in e + e~ 
annihilation was calculated in the franework of constituent quark model. The results were obtained 
for the case of final particles in pseudoscalar and vector states close to the threshold. Following the 
same procedure, we examinate the processes of e + e~ — > dd and qq — > dd, where D, are diquark and 
antidiquark, respectively (in our calculations we neglect masses of annihilating particles). Differential 
and total cross sections for the exclusive production of heavy diquark pairs in axial-axial, axial-scalar 
and scalar-scalar states are determined for the diquark composed of different heavy quarks. We 
also consider the case of axial diquark composed of equivalent quarks. We do not concern for the 
annihilation into the pseudoscalar and vector diquarks, since their production does not contribute in 
the leading order of 1/m expansion. If one supposes the fragmentation of diquarks into the doubly 
heavy baryons, the formulae derived may be useful in the calculation of cross sections for the pair 
production of baryons. 

In section 2, we describe basic points of the constituent quark model. Section 3 and 4 contains 
matrix elements, differential and total cross sections for the processes of e + e~ and qq annihilation 
into the pairs of diquarks. In section 5 numerical results are given. In Conclusion we summarize the 
consideration. 
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2 Basic points of the model 



In this paper we use the constituent quark model ||, which considers the quark masses and leptonic 
constants as the only input parameters. 

According to the model the diquark d = (Q1Q2) mass is equal to 

M = mi + ^2- 

Four-momenta of the quarks entering the diquark are given by 

1 7712 D 

k Q2 = Jf p -^ 

where P is the meson momentum. 

To represent The Fock state of diquark with two different quarks in the model, we use the principle 
of superposition for the wave packages with distributions Therefore, we find the following: 

in the case of scalar diquark 



and for the axial diquark 

where C is a matrix of charge conjugation, i and j are the colour indices, e = e^j^, e M is a polarisation 
vector of axial diquark, a and b denote the operators of quark creation, so that the axial and scalar 
diquarks in the Fock space are normilized in the following way: 

(Sl(P)\Si(P')} = (2tt) 3 M(P-P'), (3) 
(4(P,A)|4(^,A')> = {2'K)% j 6^6{P - P'), (4) 

A and A' are diquark polarisation indices. 

If the heavy quarks are identical, then the Pauli principle must be taken into account. So, the above 
formulae have to be divided by a/2 and antisymmetrized over the permutaions of a and b operators. 

The heavy quark propagator has the form 

S(k) = (k^ + m)D(k), 

where 

D~\k) =k 2 -m 2 . 
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3 Amplitudes and cross sections for the e + e annihilation 



In this section we consider the e + e~ annihilation into the diquark and antidiquark pairs. The diagrams, 
which contribute into the processes in the leading order, are shown in Figs, la and lb. The colour 
factor for the processes involved is equal to 

Colour^- = -j^ij , 

where i, j are the colours indices of diquark and antidiquark, correspondingly. 



3.1 Annihilation into pairs of scalar diquarks 

The matrix element for the pair production of scalar diquarks may be writen as follows: 



647T 2 f 

M ss = -i— J f Sy |^(0)| 2 (P; - P,% (5) 



where f ss is equal to 



fss = M {<W- S )-^ + <W^-^ 2 ) a - {S) 



/ m\ \ q 2 m 2 / m 2 n q^i \ t , ta . 



and denotes the leptonic vector current, ^ s (0) is a wave function of scalar diquark at the origin. 
The charges of the quarks Q\ and Q 2 are equal to q\ and q 2 . P', P are the four-momenta of scalar 
diquark and antidiquark, correspondingly. 

After simple algebraic calculations for the differential cross section da ss / dcos9 we get the following 
expression: 



64^^ |M> S (0)| 4 1 (1-cos 2 ^), (7) 



dcosO 3s 3 

where 6 is the angle between the momenta of lepton and diquark. 

The total cross section for the exclusive production of heavy scalar diquark pairs in e + e~ annihi- 
lation is equal to 



f 2 ( 4M 2 \ 3/2 

o- M = 2567r 3 |* s (0)| 4 (l- — ] . (8) 



3.2 Annihilation into axial and scalar diquarks 

The matrix element for the pair production of axial antidiquark and scalar diquark may be writen as 
follows: 

Mas = ^Sijfas^iO^a^e^eaPpq^l^ (9) 

where f as is equal to 



^a(O) is a wave function of axial diquark at the origin, q = P + P'. 

Carring out algebraic calculations for the differential cross section da as /dcos9 we find the following 
expression: 

da f 2 ( 4M 2x3/2 

aa as CA * J aS \,t, / n W, fn\\2 I i 



.cos* 64^0|* s (O)*„(O)|^l-_ J (2-sin^). (10) 

The total cross section for the exclusive production of heavy scalar diquark and axial antidiquark 
pairs in e + e~ annihilation is given by the following formula: 

f 2 ( 4M 2 \ 3/2 

a as = 5127T 3 |M/ s (0)M/ a (0)| 2 ( 1 - — J . (11) 

In the above formulae the difference between the masses of axial and scalar diquarks is neglected. 
3.3 Annihilation into pairs of axial diquarks 

The matrix element for the exclusive pair production of two axial diquarks may be represented as 
follows: 

M aa = -i^S tJ |* a (0)| 2 (/W(P; - P,)(e'*e) + ifJ((e'*g)e M - («,)<£)) l„ (12) 
where /W and fj$ are equal to 

For the differential cross section da aa /d cos 6* we calculate the following expression: 

dcr aa 512tt 3 |T ^ 4M 2 ^ 3/2 



* a (0)| 4 1 (A-Bcos 2 9), 



dcos6 3s 5 
where A and B have the form 

A = (/W) 2 (8 + (,7 - 2) 2 ) - - 2) + {fg)*tf + 2^), 

B = (/fl) 2 (8 + (i? " 2) 2 ) - - 2) + (/g) V - 2„), 

and 77 = s/M 2 . 

The total cross section for the exclusive production of heavy axial diquarks pairs in e + e~ annihi- 
lation is given by 

1024tt 3 , , slA ( 4M 2 \ 3/2 / t „ s 

° aa = ^- |*a(0)| 4 — J (3A-B). (15) 



s 



Actually, the behaviour of this cross section at large s is a aa ~ 1/s 3 . 
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3.4 The diquark containing two equivalent particles 

If the diquarks are composed of identical particles, the above formulae must be changed. Evidently, we 
have to take into account the production of axial diquark, only. Since the smallest angular momentum 
for the scalar diquark is equal to unit, it does not contribute in the leading order of 1/m expansion. 
In contrast, the angular momentum of axial diquark can be equal to zero. 

All formulae writen down for the annihilation into two axial diquarks remains valid except two 
form factors /W, /j^j. They must be transformed in the following way: 

/W = 2qa s (M 2 )Ma em (s), (16) 
/£] = 4qa s (M 2 )Ma em (s), (17) 

where q is the quark charge. 



4 Amplitudes and cross sections for the qq annihilation 

In this section we consider the qq annihilation into the diquark and antidiquark pairs. The diagrams, 
which contribute in the leading order, are shown in Figs, la, lb, 2, 3a and 3b. The diagram shown in 
Fig. 2 results in zero in our approach, when the quarks are on the mass shells. 
The colour factor for the diagrams shown in Figs, la, lb is equal to 

C ° l0Ur fo)(/m) = g^L> 

where m, I are the colour indices of annihilating quark and antiquark, correspondingly. It is easy to 
notice that this is a colour octet state. The colour factor for the diagrams shown in Figs. 3a, 3b is 
equal to 



C0l0Ur[!J. )(Jm) - —t^t^ - -SijSim, 

that appears to be a mixture of octet and singlet colour states. 
4.1 Annihilation into pairs of scalar diquarks 

The matrix element for the exclusive pair production of two scalar diquarks may be written as follows: 
Msg = 3|5 WQ)|2 ^f3±k P > _ jm _ (Pl^zllp^j Z|| , (18) 

where and /J^ are equal to 



5.„.„ 1 



fS = M 



+ 



mi as (s) - — [<jp s ) • ^ + <jp s ) ■ m j) as{sh (l9) 



V 



® = ^ a '^<^' (20) 

and p is the four-momentum of annihilating quark. 
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For the differential cross section da ss /d cos 9 we get the following expression: 



d<j ss 8tt 
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l*.<0)l< 2/1? + M -^ose) + <#!(. - ^)co^ ■ 



dcostf 81s 3 1 AV yi I \ Jss A JSS V s / 2 v s 

/ 4M 2 \ 3/2 

(1-cos 2 ^). (21) 



The total cross section for the exclusive production of two heavy scalar diquarks in the qq annihilation 
is given by the formula 

8tt 3 |T / 4M 2 \ 3/2 11 A 4M r 



4.2 Annihilation into axial and scalar diquarks 

The matrix element for the pair production of axial antidiquark and scalar diquark may be represented 
in the following way: 

M as = ^*:(0)* a (0) (-fS^ 'rtrPpea* - fS - 2J f^) Z„, (23) 



where /j^and ff} are equal to 



j« . M^JM.^my (s) , (24) 



-M 2 / W 2 / mfm 2 
The differential cross section da as /d cos 6* are given by 



da a , 64tt 3 



dcos# 81s 4 



1^(0)^(0)1^1- (l- 1 ^) (1 + cos 2 ^) (26) 

The total cross section for the exclusive production of heavy axial antidiquark and scalar diquark pairs 
in the qq annihilation equals 



Here we have also neglected the difference between masses of axial and scalar diquarks. 
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4.3 Annihilation into pairs of axial diquarks 

The amplitude for the exclusive pair production of two axial diquarks in the process may be represented 
as follows: 

= 32^ Q ,(4mpo M5 | / ( + } 

+ (-P'e) • ((Pp)(leT) + {lP)(pe'*)) + (ee'*) ■ ((IP)( P P') + (IP')( P P)) 
+^ • ((pe)(Je*) + (^)(pe'*)) ) • (^C - ^Mim) (28) 
+ (f£{(el)(Per) - (P'e)(e'*l)) + /W^K*)) • feC)}- 
The differential cross section da aa /d cos 9 is given by the following expression: 



• (a 4 cos 4 + 03 cos 3 + a 2 cos 2 + ai cos + a ) . (29) 

The total cross section for the exclusive production of heavy axial diquarks pairs in the qq annihilation 
has the form 



200tt 3 ^ 9x M 10 |T . ,., / 4M 2 /2 2 \ 
<r M = — — a 2 (4m 2 )a 2 (4m 2 ) 6 6 7 |fr a (0)|Vl -a 4 + -a 2 + 2a , (30) 

oi T7X^T/2 2 S V ^ \0 O / 

where we have used the following definitions: 

<Xs(s) 8mfro| ( ' a 8 {lml)m 2 a s {Aml)mi \ 

Jaa " a s (Amj)a s (Am 2 2 ) 5M 2 ^ mf m| J ' 1 j 



f[2] ^sVy , "H^'^J (on\ 



a s (s) Am\m\ ( a s (Am\) a s (4m 2 ) 
a s (4m 2 )a: s (4m 2 ) 5M \ mf ml 

a 4 = -^- S 2 ( S -4M 2 ) 2 , (33) 



400 

T 2\ 3/2 



a 3 
a 2 



1 e 3 / A/f 2 \ ' 

^(2& + /I 1 J(6M 2 + S ))^l-4^-j , (34) 

^^(s " 4M 2 ){4/W/?J S ( S - 2M 2 ) + (/W) 2 (12M 4 - 4M 2 S + S 2 ) 

+s( §(12M 6 - M\s) + 4(/[ 2 J) 2 ( S - 2M 2 ) ) }, (35) 

„2 / T\yT2 

" 1 - 4— ( 2/[ 2 J S ( S + 4M 2 ) + /W(-24M 4 + 2M 2 S + S 2 ) ), (36) 



8M 2 . 

a = --^ fej/^SM 4 - 6M 2 S + S 2 ) + (/W) 2 (-48M 6 + 28M 4 S - 8MV + s 



16M 4 

3; 

25 



+2s( ||m 4 s(4M 2 + s) + /1 2 J(-16M 4 - 4M 2 s + 2s 2 ) )Y (37) 
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4.4 The diquark with two identical particles 

In this section we consider the qq annigilation into pair of diquarks composed of two identical heavy 
quarks. So, as in e + e~ annihilation, we have to take into account the axial diquark production, only. 
The amplitude has the form 

M aa = 51 f^3 2) l^(0)| 2 {( (-PO ■ ((pP')(le) + (lP')(pe)) 

+(-P'e) ■ ((Pp)(le'*) + (lP)(pe'*)) + (ee'*) ■ ((lP)(pP') + (lP')(pP)) 
+ | " ((pe)(le'*) + (!e)(pe»)) ) • (t?tf m - ^Mim) 

+ (fS{(el)(Per) - (P'e)(e'*l)) + fS}(Pl)(etf)) ■ fee)}- (38) 
The differential cross section da aa /d cos 6 equals 



■P^-a = 211 2 !f ilfrV aW \\ 1 - —(a* cos 4 9 + a 3 cos 3 9 + a 2 cos 2 + a, cos + a ). (39) 
ct cosy 81 M z s' V s 

The total cross section for the exclusive production is given by 

= 2 -grik^ 1 ^ 1 V 1 " T (s 04 + 3° 2 + 2a °) ' (40) 

where we have used the following definitions: 

m = M2 (4i) 

Jaa " a s (M 2 ) 5 ' y ! 

m = (42 ) 

Jaa a s {M 2 ) 5 ' V ' 



5 Numerical estimates 

The cross section ratios of the exclusive heavy diquark pair production to the respective heavy quark 
cross section in e + e~ annihilation, 

a(e+e- -> QQ) = = J 1 - 1 + , 

are shown in Figs. 4a, 4b. 

We see that the most optimistic expectations for the pair production is given for the cc-diquarks at 
B-factories, when the hogh luminosities can yield several thousands pairs of doubly charmed baryons. 

The cross section ratios of the exclusive heavy diquark pair production to the respective heavy 
quark cross section in qq annihilation, 

8na 2 F 4mX/ 2m 2 \ 
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are shown in Figs. 4c, 4d. 

Figs. 4a, 4c represent the case of diquark composed of different quarks and Figs. 4b, 4d give the 
case of diquark composed of identical quarks. We put 

m c = 1.55 GeV, 

m b = 4.9 GeV, (43) 
A QCD = 0.2 GeV. 

The values of diquark wave functions at the origin have been taken from [1[ • 

In the quark-antiquark annihilation the pair production of doubly charmed baryons gives about 
10~ 5 fraction of charm yield at ^Js < 100 GeV in the hadron collisions, while at higher energies the 
gluon-gluon fusion will dominate. In experoments at fixed target the threshold behaviour results in an 
additional suppression. 

The inclusive production of doubly heavy baryons was considered in [|J, so that we see that the 
pair production results in the suppression factor about 0.1. 



6 Conclusion 

In this paper we have considered the exclusive production of doubly heavy diquark pairs for the axial- 
axial, scalar-scalar and axial-scalar states in the framework of the constituent quark model. The 
matrix elements, differential and total cross sections for the processes of e + e~ and qq annihilation have 
been given. We have calculated the pair production of diquarks composed of identical heavy quarks. 
The obtained formulae can be used in the calculation of cross sections for the diquark fragmentation 
into the doubly heavy baryons in the processes of e + e~ annihilation and proton- antiproton inelastic 
collisions. 

We have found that the yield of doubly heavy baryon pairs could really reach 10 3 -hl0 6 depending on 
the energies and luminosities of accelerators, e + e~ colliders and fixed target experiments with hadron 
beams. 

The authors thank Dr. V.V.Kiselev, who asked us for the calculations of pair production of doubly 
heavy diquarks and explained some points in the problem. 

This work is in part supported by the Russian Foundation for Basic Research, grants 99-02-16558 
and 00-15-96645. 
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Figure Captions 

Fig. 1. The diagrams of e + e~ (qq) annihilation into the pair of doubly heavy di quarks due to the single 
photon (gluon) exchange. 

Fig. 2. The qq annihilation into the diquarks due to the three-gluon interaction. 

Fig. 3. The diagrams of double gluon annihilation of light quarks into the pair of doubly heavy diquarks. 

Fig. 4a. The ratios of total cross sections for the exclusive production of be diquark pair and the cc 
production in e + e~ annihilation for the axial-axial, axial-scalar and scalar scalar states. 

Fig. 4b. The ratios of cross sections for the exclusive production of bb and cc diquark pair and the 
respective bb and cc production in e + e _ annihilation. 

Fig. 4c. The ratios of cross sections for the exclusive production of be diquark pair and the cc production 
in qq annihilation. The notations are the same as in Fig. 4a. 

Fig. 4d. The ratios of cross sections for the the exclusive production of bb and cc diquark pair and the 
respective bb and cc production in qq annihilation. 



11 



Fig. la 



Fig. lb 
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